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ABSTRACT 

Strain-displacement  relations  for  thin  shells  ralid  for  large 
displacements  are  derived.  With  these  as  a  starting  point  approximate 
strain-displacement  relations  and  equilibrium  equations  are  derived  by 
making  certain  simplifying  assumptions.  In  particular  the  middle  surface 
strains  are  assumed  small  and  the  rotations  are  assumed  moderately  small. 
The  resulting  equations  are  suitable  as  a  starting  point  for  stability 
investigations  or  other  problems  in  which  the  effects  of  deformation  on 
equilibrium  cannot  be  ignored,  but  in  which  the  rotations  are  not  too 
large. 

The  linearized  forms  of  several  of  the  sets  of  equations  derived 
herein  coinelde  with  small  deflection  theories  in  the  literature. 
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INTRODUCTION 

One  of  the  Important  uses  to  which  a  large  displacement  theory  of 
thin  shells  can  be  put  is  the  investigation  of  stability.  The  many  papers 
in  the  literature  on  the  stability  of  shells  have  dealt  almost  exclusively 
with  cylinders,  spheres,  and  cones,  and  the  differential  equations  governing 
the  phenomenon  have  been  derived  specifically  for  these  geometrical  shapes. 
It  would  seem  to  be  desirable  to  have  a  unified  treatment  based  on  a  general 
theory  for  an  arbitrary  middle  surface.  The  practically  important  cases  of 
the  shallow  shell  and  the  shell  of  revolution  with  symmetric  deformations 
have  been  adequately  treated  (references  1  and  2)  but  the  general  problem 
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presents  difficulties  not  found  in  the  special  cases.  It  is  the  purpose  of 
the  present  paper  to  derive  an  exact  theory  for  large  deflections  of  a  thin 
shell  with  an  arbitrary  middle  surface  and  then  by  making  certain  simplifying 
assumptions  to  derive  from  this  several  theories  suitable  for  application. 

An  incomplete  treatment  of  the  general  large  deflection  theory  of  thin 
shells  has  been  given  by  Novoahilov  in  reference  3.  He  derives  a  theory  for 
small  middle  surface  strains  but  does  not  go  into  detail  on  further 
simplifications  or  discuss  approximate  squilibriua  equations.  He  indicates 
that  the  next  step  is  to  assume  that  the  rotations  are  small.  His  results 
for  the  small  strain  theory  differ  from  those  in  the  present  paper. 

In  the  western  literature  there  have  been  several  papers  dealing  with 
the  general  large  deflection  theory,  in  particular  the  paper  by  Synge  and 
Chien  (reference  li),  the  series  of  papers  by  Chien  (references  5  and  6)  and 
the  recent  paper  by  Ericksen  and  Truesdell  (reference  7).  The  intrinsic 
theory  of  shells  developed  by  Synge  and  Chien  avoids  the  use  of  dlsplace- 
msnts  as  unknowns  of  the  problem.  The  theory  of  shells  is  deduced  from  the 
three  dimensional  theory  of  elasticity  and  then  by  means  of  series 
expansions  in  powers  of  a  small  thickness  parameter  approximate  theories 
of  thin  shells  are  derived.  A  large  number  of  problem  types  is  found 
classified  according  to  the  relative  orders  of  magnitude  of  various 
quantities.  This  approach  has  been  criticized  by  severtl  authors  (see 
references  8,  9  and  10),  In  the  paper  of  Ericksen  and  Truesdell  there  is 
a  unified  treatment  of  shells  and  curved  rods  developed  as  two  and  one 
dimensional  theories  respectively  without  an  attempt  to  deduce  these  from 
the  three  dimensional  theory  of  elasticity.  The  constitutive  relations 
ere  purposely  left  out  of  consideration  since  they  are  unnecessary  for 
the  description  of  strain  and  the  establishment  of  equilibrium  conditions. 

The  authors  are  interested  only  in  exact  theory  and  do  not  discuss  the 
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simplifications  of  small  strains  and  rotations. 

The  papers  of  Synge  and  Chien  and  that  of  Ericksen  and  Truesdell  on  the 
shall  problem  are  unorthodox  and  difficult  to  relate  to  most  of  the  rest  of 

i 

the  literature.  The  authors  have  aimed  at  a  maximum  of  generality,  perhaps 
more  than  necessary  for  the  technological  applications  of  the  theory. 

In  the  present  paper  a  large  deflection  theory  for  thin  shells  is 
developed  in  which  transverse  shear  and  normal  strains  are  neglected.  With¬ 
out  further  approximations  the  equations  are  very  complicated,  but  sine* 
large  middle  surface  strains  are  almost  never  encountered  in  engineering 
applications,  an  effort  is  made  to  derive  simpler  equations  based  on  the 
assumption  of  small  strains.  Considerably  more  simplification  is  gained 
when  the  rotations  as  well  as  the  middle  surface  strains  are  assumed  to  be 
small.  In  this  way  a  theory  for  shells  is  derived  comparable  to  the 
von  Karman  theory  for  plates.  The  equations  corresponding  to  other  minor 
simplifications  are  also  presented.  The  main  part  of  the  paper  is  written  ■ 
in  tensor  notation  but  the  principal  results  are  reproduced  in  the  ordinary 
notation  in  an  appendix. 


DERIVATION  OF  EXACT  EQUATIONS  FOR  LARGE  DISPLACEMENTS 
Geometrical  Preliminaries 

Let  the  undeformed  middle  surface  of  the  shell  be  given  by  the  equations 
x1  -  x\za)  (i  -  1,  2,  3  }  a  -  1,  2)  (1) 

i  JCL 

where  the  x  are  cartesian  coordinates  in  space  and  the  are  curvi¬ 
linear  coordinates  on  the  surface.  Let  the  displacements  of  material 

points  on  the  middle  surface  of  the  shell  be  resolved  into  components 
tangential  and  normal  to  the  undeformed  middle  surface  as  expressed  by  the 
following  equation 


U1  -  u“(5)xi„  +  W(Z)Nl  (2) 

vher«  xd  ■  ——  are  tangent  v-3-jtora  to  the  coordinate  curves  on  the  un- 

'*  3«“ 

deformed  middle  surface  and  H*  is  the  unit  normal  to  the  undeformed 
middle  surface.  In  this  paper  the  coordinates  5®  will  be  used  to  label 
material  particles  on  both  the  undeformed  middle  surface  and  the  deformed 
middle  surface. 

Some  of  the  important  formulas  of  the  theory  of  surfaces  vill  be  used 
repeatedly  and  are  reproduced  here  for  convenient  reference.  For  the  un¬ 
deformed  middle  surface  the  formula  for  the  squared  element  of  arc  ds2 
in  terms  of  the  first  fundamental  form  g  .  is 

°p 

ds2  *  x±,aX^  '  gcf3  <3> 


the  element  of  area  is 


dA  «Vg  d^d«2 


where  g  is  the  determinant  of  g^  j  and  the  equations  of  Oauss, 
Waingarten  and  Codas*!  are 


N1  -  bjx1 
jO  a  ,p 


bap,Y  '  V.P  ^ 

where  a  comma  denotes  covariant  differentiation  with  respect  to  the  metric 
g^  and  where  the  second  fundamental  form  b^  ,  as  here  used,  differs 
in  sign  from  the  usual  definition. 

After  the  displacement  given  by  equation  (2),  the  material 
particle  originally  at  xd  will  move  to  the  point  y*  given  by 
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y1  =■  x1  +  O1 

.  x1  +  U“x-  +  KN1 
»« 

(8) 

» 

• 

• 
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This  is  the  equation  of  tha  deformed  middle  surface  in  terms  of  the 

parameters  .  Tangent  vectors  to  the  coordinate  curves  on  the  deformed 

middle  surface  are  given  by 

..  ■  i 

>;  •'  » 

•.  | 
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V  -  XeKr  + 

(9) 

'  \ 

1 

where 

• 

Xof?  "  K/,?  “  4*  +  V  +  V 

(10) 

*  - » 

«• 

(11) 

V 

e 

ti  -  y1  n1  .»  .A, 

a  ,a  |(x  c  p 

(12) 

i » 

:  l 

Also  define 

i‘  : 

v  -  x1/1 

(13) 

•  » 

• 

cos®  ■ 

(11*) 

• 

where  is  the  unit  normal  to  the  deformed  middle  surface  (see  equation  - 

(73)).  Indices  on  X^  ,  na  and  va  will  always  be  raised  or  lowered  with 

the  metric  g.  .  Where  necessary  a  bar  indicates  a  quantity  defined  with 
°P 

’  V 

• 

respect  to  the  deformed  middle  surface.  The  squared  element  of  arc  on  the 

e 

deformed  middle  surface  is  given  by 

• 

ds2  -  y^p  -  0^  dC“d«P 

In  terms  of  X_«  and  u 

up  a 

(15) 

• 

• 

°ep  *  XaXrP  +  Wf> 

The  element  of  area  dX  on  the  deformed  middle  surface  is  given  by 

(16) 

• 

• 

• 

- - - - -  - - -  — -  - 

• 

dK  ■-/3  d^d£^  "f^g  8A 

The  equations  of  Oauss,  Weingarten  and  Codazsi  are 

(17) 

• 

(18) 

*.  •  **. 

(19) 

« 

Bap/r  “  V/P 

(20) 

where  B.  is  the  second  fundamental  form  of  the  deformed  middle  surface. 

“P 

X  slash  is  used  to  denote  covariant  differentiation  with  respect  to  the 
metric  0  „  .  For  scalars  and  cartesian  tensors  (as  W  ,  y*  or  )  there 

“p 

is  no  difference  between  a  slash  and  a  comma,  in  expression  for  B  .  in 

“P 

terms  of  X  „  ,  u  ,  v  and  cos«>  can  be  derived  as  follows 

Op  3  u 

v  ■  *lA  -  ■  W 

■  *  Kr'lW’1’1 

m  (b?V^,p)c<,8“  ”  (Xra,P+V«)yY  (a) 

The  following  identity  will  prove  to  be  useful 

yja®1  ■  »%<,  +  \  co*“  -  0  (22) 


In  the  coordinate  system  of  the  deformed  middle  surface  the  equilibrium 
equations  of  the  shell  are  the  same  as  in  the  linear  theory  and  need  not  be 
derived  here.  They  are  (see  reference  10) » 
force  equilibrium 


+  B^Q°  +  pP  ■  0 

(23) 

-  *  P  "  0 

(2U) 

). 

- 1  - 

moment  equilibrium 

• 

i- 

1 

o 

« 

°?y 

(25) 

I. 

4 

t  Jr* 

\ 

ii 

s-  (N°p+B yP)  -  0 

o P  Y 

(26) 

'.M 

■f 

) 

where  N8^  is  the  membrane  stress  resultant,  M®P 

is  the  bending  moment 

•s* 

0 

»  I" 


resultant,  and  Q°  is  the  transverse  shear  stress  resultant,  all  defined 
with  respect  to  the  deformed  shell.  The  quantities  p®  and  p  are  applied 
load  intensities  per  unit  of  area  of  the  deformed  middle  surface,  e_D  is 

op 

the  covariant  permutation  tensor  in  the  deformed  coordinate  system. 

The  above  equilibrium  equations  are  exact  tut,  of  course,  the  ten 
stress  quantities  entering  into  them  do  not  furnish  a  complete  description 
of  the  state  of  stress  throughout  the  thickness  of  the  shell.  However,  in 
thin  shell  theory  it  is  always  assumed  that  the  state  of  stress  is  adequately 
described  in  terms  of  these  quantities. 

?inlte  Strains 

The  strain  quantities  entering  into  a  thin  shell  theory  are  a  matter 
for  definition.  The  literature  of  the  subject  shows  a  wide  variety  of 
ohoices  of  strain-displacement  relations,  particularly  for  the  bending 
strains.  Some  choices  have  been  shown  to  be  better  than  others  (see 
reference  11)  but  at  the  present  time  no  set  of  conditions  sufficient  to 
render  the  ohoice  unique  has  been  generally  agreed  upon.  In  the  present 
piper  the  ohoice  has  been  guided  by  two  considerations,  the  first  of  which 
was  the  desire  to  derive  a  theory  which  admits  a  principle  of  virtual  work. 
This  requirement  forces  a  close  relation  between  the  equilibrium  equations 
and  the  strain-displacement  relations.  The  second  consideration  was 
simplicity.  The  resultant  cnoice  will  be  shovi.  to  furnish  an  adequate 
description  of  the  deformation  of  the  shell  provided  the  Kirchoff  hypotheses 
are  accepted  as  adequate  descriptions  of  the  displacements. 
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Let  S  be  a  simply  connected  region  on  the  deformed  middle  surface 
enclosed  by  the  curve  J  .  The  following  identity  follows  from  equations 
(23)  to  (26). 

J  ♦  (Q/a-B^H^+p)65 

S 

+  (M^-Qp)60p  +  S^(N^+B“M^)6^dX  -  0  (27) 

By  application  of  the  divergent'  theorem  for  a  curved  surface  (27)  may  be 
transformed  into  the  following  identity  which  is  the  preliminary  fom  of  the 
principle  of  virtual  work  and  all  subsequent  derivations  will  proceed  from  it. 

J>  (N^fiUp  +Qa6V7+Map6jZfp)nt  ds  +  J  (p^fiffp+pSlDdl 

Z  S 

5 


+  MaP(«0p/c-JYpB^)]dA  (28) 

where  the  virtual  displacements  6^  ,  and  rotations  6 0a  ,  50  refer 
to  components  in  the  directions  of  the  tangents  and  normal  to  the  deformed 
middle  surface.  In  (28)  the  terns  on  the  left  hand  side  are  interpretable 
as  the  external  virtual  work  of  edge  loads  and  surface  loads  respectively. 
The  right  hand  aide  of  (28)  might  be  interpreted  as  internal  virtual  work 
if  the  coefficients  of  N°^  ,  Q°  and  were  identified  with  strain 
increments.  Such  an  identification  will  be  postponed  until  later.  First 
these  coefficients  will  be  written  in  a  different  form. 

By  definition 

6U1  -  SU^y1  +  6ft  jf1  -  6Uax^\  +  6W  N1  (29) 

Now 
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(30) 


6U)P  -  (eW/p-B^sUV  ♦  (6^+Bj6VJ)yJo  -  6U*p 

■  (6w,p-v6D<I)Nl  +  (60°p+bp6W>xi« 

From  which: 


5V« +  •  vx* +  v*«  (31) 

&/a-  *  vy6x£  *  cow  6na  (32)' 

The  rotation  around  the  nomal  60  is  given  in  terns  of  dieplaoemente  by 
the  fonmila 

60  -  |  e^aUp/a  (33) 

By  the  use  of  (31)  and  the  fact  that  B  Q  is  symmetric  (33)  becomes 

®P 

60  -  |  eaP(Xrp6X^p6no)  (31*) 

Also  60  -  ^VX^V^X^XJ-H^)  (35) 

From  (31)  aid  (35)  the  coefficient  of  N0^  in  (28)  is 

i<Vl>pWh8V  ’  1 6(V>«V  - 1  60of*  (36) 

The  natural  definition  of  the  finite  membrane  strain  is  thus 

*«*  -  t<  W  (37) 

The  coefficient  of  Qa  in  (28)  is 

K  ‘  6ff/a  -  4®^  +  %  "  \6Xl  *  C0Sa  +  {38) 

Since  the  intention  is  to  derive  a  theory  in  which  transverse  shear  strains 

are  neglected,  set  6y  »  0  which  gives 
ct 

60a  *  -  v  «xj  -  costa  6(*a  (3?) 


which  serves  to  relate  rotations  to  displacements.  From  (9),  (13),  (ll*)  and 
(39) 
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5y'  -  ... 

J  yd  Y  a 

or  since  ffV1  ■  0  it  follows  that 


ifoy1  -  v  8X^  +  cofu  6u  »  -  60 


60a  -  y^ff1 


m 


(la) 


A  finite  transverse  shear  strain  ra  consistent  with  (38)  and  the  require¬ 
ment  ■  0  nay  be  defined  as  follows 


*  V  v  +  a  cos<a 
ay  a 


(U2) 


The  coefficient  of  m“P  in  (28)  may  be  found  in  terms  of  B  and  G 

ap  q 

as  follows.  From  (Id) 


■  ^ A 


(U3) 


sinoe  M-  -  0  ,  Now  recall 


**  ■ 


so 


Which  gives 


From  (9)  and  (35) 


“op  ■  •  «>* 

■  «v.  *  *yy,t 

»v  ■  6v  • 


-  Kss'»v 


(WO  ‘ 


(1*5) 


From  (U5)  and  (U6)  the  coefficient  of  in  (28)  reduces  to 
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?•  •  < 


(1*6) 
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SB. 


BY6E„ 

a  0y 


(U7) 


Using  the  foregoing  results  the  right  hand  side  of  (28)  nay  be  written 

f  [H^SEaf  *  (k8) 

S 


There  is  obviously  some  difficulty  in  defining  a  finite  bending  strain 
tensor  because  the  coefficient  of  in  this  expression  is  not  the 

exact  variation  of  anything.  However,  a  way  to  proceed  suggests  itself  if 
(lj8)  is  rewritten  in  the  following  form 

J  [(H^-B^SE  +  M^SB^ldA  (U9) 

2 


Define  a  finite  bending  strain  tensor  by 

Kc$  "  V  "  V 

Define  a  modified  membrane  stress  tensor  by 


(50) 


-  N°P  -  b^my°  (51) 

Sinoe  B  .  ie  symmetric  there  will  be  no  loss  in  generality  by  defining  a 
°P 

modified  bending  moment  tensor  by 


a*  -  |(MaP+rfa)  (52) 

Note  that  the  third  moment  equilibrium  equation  (equation  (26))  is 
equivalent  to  the  statement  that  is  symmetric.  In  terms  of  the 

newly  defined  quantities  ( U9 )  becomes 

j  or^A^di 

s 


(53) 


The  details  will  not  be  shown  here  but  this  expression  for  the  internal 
virtual  work  aay  be  derived  from  the  three  dimensional  theory  by  integration 
through  the  thickness  of  the  shell  and  without  approximation  provided  the 
displaeanents  are  restricted  by  the  Kirchoff  hypotheses.  The  appearance  of 
the  modified  tensors  V*  and  rather  than  K°^  and  may  seem 

somewhat  strange  but  then  there  is  no  reason  why  these  quantities  should  not 
be  adopted  as  the  stress  quantities  entering  into  the  theory  rather  than 
and  In  fact  there  is  an  advantage  in  that  both  and 

are  symmetric  and  thus  there  are  fewer  unknowns  to  deal  with.  This  must 
simply  mean  that  the  equilibrium  equations  containing  and  M8^  are 

slightly  more  general  than  is  appropriate  for  a  theory  in  which  the  dis¬ 
placements  are  restricted  by  the  Kirchoff  hypotheses.  Moreover,  in  the 
present  theory,  there  are  the  same  number  of  stress  quantities  as  there  are 
strain  quantities.  If  both  a  principle  of  minimum  potential  energy  and  a 
principle  of  minimum  complementary  energy  (for  an  elastic  shell)  are  to  be 
possible  of  formulation  in  the  theory,  then  it  is  ordinarily  necessary  that 
the  constitutive  relations  be  invertible.  This  is  possible  only  if  there 
are  the  same  number  of  stress  quantities  as  strain  quantities. 

The  f&o  strain  tensors  E  .  and  K  .  (or  equivalently  the  two  tensors 

up 

0^  and  )  furnish  an  adequate  description  of  the  deformation  of  the 
■hell  as  shown  by  the  following  argument.  In  the  first  place  the  deformation 
of  the  shell  is  completely  described  in  terns  of  the  displacementsof  points 
on  the  middle  surface  provided  the  displacements  throughout  ths  thickness  > 
are  restricted  by  the  Kirchoff  hypotheses.  Secondly,  from  the  theory  of 
surfaces,  a  knowledge  of  0^  and  as  functions  of  S8  and  subject 

to  the  Oauss  and  Codazzl  integrability  conditions  v — ch  in  the  present  case 
are  equivalent  to  compatibility  conditions)  completely  determines  the  deformed 
middle  surface  together  with  a  coordinate  system  (the  deformed  system) 


'V 

A 


'V# 


except  for  a  rigid  body  potion. 
Modified  Equilibrium  Equations 


Since  new  stress  quantities  have  been  introduced  the  equilibrium 
equations  (23)  to  (26)  are  no  longer  appropriate  for  the  theory  being 
developed  here.  Appropriate  equations  can  be  derived  from  the  expression 
(53)  for  internal  virtual  work. 

J  ^v^Vd*  *  J  n*(8V«4V0) 

5  J 

+  ^(-6%^/86irY+2Bj6irY/a+B^6W)]dr 

-  f  *  ifgfiV?  ♦  ds 

C 

-Jt^B^+BP/ar«)6Up 

J 

♦  <&> 

The  line  integral  around  C  is  the  external  virtual  work  of  the  edge 
forces  ani  moments.  If  a  principle  of  virtual  work  is  required  to  hold, 
and  if  tH  portion  of  the  shell  within  7  is  in  equilibrium,  then  the 
internal  virtual  work  must  equal  the  external  virtual  work  for  arbitrary 
virtual  displacements.  Thus  the  condition  of  equilibrium  is  that  the  last 
integral  in  (5U)  vanishes.  This  leads  to  tho  following  equilibrium  equations 
(surface  forces  have  been  omitted  in  this  derivation  for  simplicity). 

+ + -  °  (») 

%  -  v*  -  Biy^  ■  °  (56) 

Equations  (55)  and  (56)  are  in  fact  identical  to  the  equations  (23)  to  (25) 
with  Qa  eliminated.  Equation  (26)  is  accounted  for  by  the  symmetry  of 


fl“P  .  If  73“  defined  by 


7f  -  ^  (57) 

is  introduced  as  an  approximation  to  Q°  ,  then  the  equilibrium  equations, 
in  an  expanded  form,  may  be  written 


o 

■ 

(58) 

73“  -  b  X*  -  b\  If#  -  o 
/«  a?  a  PY 

(59) 

O 

• 

(60) 

Tha  principle  of  virtual  work  reads 

6  ♦  Qa6\7  +  da 

J  Y  p  p  0 

C 

■  J  +?f(64+6VBa6°P) 

.  J  ((J^SE^+Tf^^fiK^dT  (6l) 

3 

It  is  something  of  a  matter  of  personal  preference,  but  it  is  this  form  of 
the  principle  of  virtual  work  (with  the  73“  and  fa  terms  present)  which 
will  be  used  in  the  remainder  of  the  paper. 

Sim  STRAIN  APPROXIMATION 

Equilibrium  Equations 

In  almost  all  practical  applications  of  shell  theoiy  the  middle  surface 
strains  are  small  whether  or  not  the  displacements  are  small,  in  which  case 
the  exact  equations  can  be  simplified  somewhat.  Since  by  definition  the 
middle  surface  strain  is 


small  middle  surface  strain  means  that  the  intrinsic  geometiy  of  the  deformed 
middle  surface  is  almost  the  same  as  that  of  the  undeformed  middle  surface. 
Thus  for  sufficiently  small  middle  surface  strains  (and  exactly  for  in- 
extensional  bending)  covariant  differentiation  with  respect  to  G^  may  be 
replaced  by  covariant  differentiation  with  respect  to  g^  .  The  equilibrium 
equations  become 

1*  *  ♦  0$^,.  ♦  ?  -  0 

-  BIV^  +  P  *  0  (62) 

VP-tfmO 

>a 


Strain-Displacement  Relations 

A  short  Investigation  shows  that  the  expression  for  E^  given  by  (37) 

and  (16)  is  not  simplified  at  all  by  the  assumption  of  small  strains.  The 

expression  for  given  by  (21)  (and  given  by  (SO))  is  slightly 

simplified  because  the  expressions  for  cosoo  and  vB  in  terms  of  X^ 

and  (i  (see  equations  (7k)  and  (7S)  later  in  the  paper)  can  be  simplified 
c 


oj  T9 plac^g  V®  ^  ^-9  expression  for  is  obviously  still 

quite  complicated. 

Constitutive  Relations 

Consistent  constitutive  relations  for  the  linear  small  strain  theory  of 
thin  elastic  shells  have  been  derived  in  references  11,  12  and  13.  These 
derivations  require  only  minor  modifications  in  the  case  of  finite  dis¬ 
placements  and  small  strains  so  they  will  not  be  reproduced  here.  ?or  a 
thin  shell  of  uniform  thickness  h  composed  of  an  Isotropic  hookean 
material,  the  constitutive  relations  are  the  same  as  in  Love's  first 
approximation,  namely  the  linear  relations 


(63) 


EhS^  -  U+v^gp/-6  -  vg^g^6 

n  *  (1+v)gaYgp6^6 "  Vr&w5 

According  to  reference  11  these  relations  may  be  used  even  if  the  definition 

of  K  .  in  terms  of  displacements  is  altered  by  the  addition  of  terms  of  the 
°P 

form  BrE  .  .  A  similar  argument  to  that  in  reference  11  shows  that  In* 

arp 

may  be  altered  by  addition  of  terms  of  the  form  B^Jr^  .  Note  that  in  the 
case  of  small  strain  the  indices  on  and  may  be  raised  and 

lowered  with  the  metric  instead  of  0^  with  negligible  error. 

If  the  material  of  the  shell  is  not  elastic  and  isotropic  the  relations 
(63)  and  (61*)  must  be  replaced  by  others  appropriate  for  the  material.  How¬ 
ever,  the  strain-displacement  relations  and  the  equilibrium  equations  given 
previously  are  unaffected  by  the  material  so  long  as  transverse  shear  and 
normal  strains  can  be  neglected. 

APPROXIMATION  OF  SHALL  STRAINS  AND  MODERATELY  SMALL  ROTATIONS 
The  exact  theory  was  considerably  simplified  by  the  assumption  that 
the  middle  surface  strains  are  small,  but  the  equations  are  still  very 
complicated.  Considerable  additional  simplification  can  be  achieved  if 
the  rotations  are  assumed  small  also.  This  simplification  will  be  carried 
out  in  the  following. 

For  infinitesimal  displacements  and  rotations  it  is  evident  from  (33) 
and  (38)  that  the  rotations  are  given  by  the  formulas 

*  •  I  <«> 

and 

K  ■  -w,« +  ❖p  *  -*«  <66> 

For  small  but  finite  rotations  it  is  convenient  to  think  of  the  expressions 


*  J. .  - 

ti  (65)  and  (66)  as  stations  (just  as  in  the  linear  theory  of  shells). 
Purely  for  convenience,  suppose  that  the  coordinates  f  have  the  units  of 
length  so  that  0 , 0g  and  are  dimensionless.  The  following  order 

of  magnitude  assumptions  will  lead  to  a  theory  one  step  beyond  the  linear 
theory  in  refinement. 

0  0r  §<  W  "  °(e)  (67) 

0a  or  n„-0(e)  (68) 

+  V  "  °{e2)  (69) 

where  e  is  a  number  small  compared  to  unity.  Write  X  .  in  the  form 

°P 

V  ’  gop  +  +  V+2V) 

‘  gap  "  V  +  (70) 

From  (16)  and  (37)  we  find  that  is  0(e2)  end  is  given  approximately 
by  the  expression 

\  '  K,p+Up,a)  +  V  +  I  Vp  +  \  (71) 

The  order  of  magnitude  assumption  (69)  was  made  so  that  these  terms  would 

not  domirfte  the  expression  for  E^  ,  otherwise  the  linear  theory  would 
result. 

In  order  to  simplify  the  expression  for  expressions  for  v  and 
cos«  are  needed.  An  expression  for  IT1  may  be  found  ty  taking  the  cross 
product  of  the  tangent  vectors  y*  . 

V*  ’  eij^p  (72) 

Chitting  the  details,  this  leads  to  the  following  expression  for  J71  in 

terms  of  \  and  p 
a 


-XU- 


”  ]/§  *<^eY5<‘|  X-YaX6pf,i+XrP,1«Xi5^ 


From  (13) ,  (lit)  and  (73) 


•*  -yl  av 

•  i  j/§  *“f*Y\«x«p 


(73) 

(7U) 

(75) 


From  0^  -  g^  -  0(e2)  it  follows  that  -y/§  »  1  ♦  0(*2)  j  then  from  (67) 
to  (70)  and  (,7k)  and  (75)  it  follows  that 

v8  -  -n8  ♦  0(e2)  (76) 

cosco  »  1  +  0(t2)  (77) 

From  the  foregoing  and  (21)  the  first  approximation  to  B  .  -  b  .  is 


Tap  "  Bap  "  bcp  S  2  bp^Y,o_Uo,Y^  “  ^a,p 


(78) 


A  difficulty  here  is  that  this  expression  is  not  symmetric  in  a  and 
p  .  However,  it  can  easily  be  shown  that  the  antisymmetric  part  is 


WK  ' 
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is  negligible  (compared  to  +  Tpa).  The  symmetric  part  of  can 
be  taken  as  the  first  approximation  to  the  bending  strain,  giving 
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negligible. 

4*’ 
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r 

hifV'M  >-W<w 

(79) 

f. 

• 

Now 

*  v  ■  0(‘!) 

% 

-y 

:/ 
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so 

i<W  ■  -v  * ol,2)  ’ 

(80) 

» 

• 

and  then 

Tap  ~  Tpa  *  (bpVbabrP)W  +  order  terns  » 

(81) 

* 

but  the  first  term  on  the  right  hand  side  of  (81)  vanishes  so 

Tap  "  Tpa 

• 

S 

i 


-x/- 


op 


or 


K  .  ha  „+0„  )  +  i  b£(U  -0  )  +  r  bT(U  .-VB  )  (82) 

op  r'a, p  yp,ay  u  p'  r>8  o,r'  «  <0  y>P  P»r 

Alternatively  it  can  do  no  harm  to  include  soma  higher  order  terms  in 

K  if  for  some  applications  it  simplifies  the  equations.  The  following 
°P 

approximation  of  may  be  used  instead  (see  equation  (21)) 

Kcp  ‘  bpXr«  "  ^a.p  “  bop 

-  0  „  +  blU  +  bXb  W 
*“»p  P  P  r« 


*  »!/,  *  1?t,P  *  bP°r,.  *  bpV” 


(83) 


which  is  symmetric  as  it  stands.  Both  alternative  expressions  are  linear 
in  the  displacements.  The  expression  (82)  for  the  bending  strain  is  the 
same  as  the  one  derived  in  references  11  and  llu  The  expression  (83)  does 
not  seem  to  have  appeared  before  in  the  literature.  Either  expression 
belongs  to  a  consistent  shell  theory  in  the  sense  of  ICoiter  (reference  11). 
When  transverse  shear  strains  are  neglected,  the  condition  8r  ■  0 

U 

serves  to  relate  60a  to  the  displacements  as  in  (39).  By  analogy  the 
following  Expression  should  be  nsed  for  y8  in  the  present  case 


K  *  \  ■  K  *  .  -  b?T 


(81») 


Approximate  Equilibrium  Equations 

Approximate  equilibrium  equations  corresponding  to  the  approximate 
strain-displacement  equations  may  be  found  by  the  same  method  used  to 
derive  the  equilibrium  equations  (53)  and  (36).  Since  the  strains  are 
small  dX  may  be  replaced  by  dA  and  we  have 
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J  l^\p  +  b«pW  +  +  2  e^eY65U6,r] 


• +  **£*«,/»  +  I  bI(6Ur,r6UP.r)]  +  *t*u«K,m-*;*y)  }  «* 

■  M  {0?.  ♦  ♦  |  «*<*{>„  *  WTi.r  -  §<*#%  ' 


*  #>p  +  Pfp  -  ♦  (^)tp]6W  ♦  [K^P  -  *“]60a  JdA  (85) 

Ijjr  Inspection  the  equilibrium  equations  are 

1*  ♦  ^  *  |  *  !<#%  -  &T\,  *  &  *  p"  •  0 

(»«) 


-  V*  *  <‘,/f>,|!  *  '  ■  0 


(87) 


S"jj  -  tf  -  0  (88) 

»p 

where  the  load  terms  pa  and  p  have  been  supplied.  Thesa  equations  express 
equilibrium  of  forces  and  moments  in  directions  parallel  to  the  tangents  and 
norual  to  the  undeformed  middle  surface.  In  the  left-hand  side  of  (85)  0a 
could  be  expressed  in  terms  of  displacements  and  the  term  3*5^  could 
be  omitted.  The  result  for  the  equilibrium  equations  would  be  (86)  end  (87) 
with  72®  eliminated  by  means  of  (88).  If  the  expression  (83)  is  used  for 
the  bending  strain,  the  first  two  equilibrium  equations  are  slightly  different 
and  read  as  follows 


*5  *  tyy  •  |  *  o^m,.  *  *  /  ■  0  <86)' 

v,  ■  .  <»„ir%  -  .  p  -  0  (67)' 

Boundary  Conditions 

The  Kirehoff  boundary  conditions  may  be  obtained  from  the  boundary 


integral  in  (85)  which  when  written  out  reads 


}  |CB°P  -  |  *  |  bV?  -  |  ♦  (lf31%a)6W  +  Tf^60a  ]  np  ds 

°  (8?) 
Let  t°  be  the  unit  tangent  to  the  curve  C  ,  then  n^  •  s^t’*’  is  the  unit 
normal  to  C  in  the  surface  S  .  Let 

K  ■  *  K\  {90) 

where  0  and  0„  are  scalars.  From  (?0) 
ra  n 


0  K  0  ta  -  (-W  +bYtf  )t“ 
,a  a  y 

.  -  ♦  bYU  t° 

ds  a  y 


(91) 


Obviously  08  is  not  independent  of  W  and  tfa  on  C  .  The  last  tern 
in  (89),  namely, 

f  A  *>  +  t£y6)ta  ♦  6ya]np  ds  (92) 

c  c 


becomes,  upon  integrating  by  parts, 

f  f^(T5“Ptanp)6W  +  0  +  ^Vp6^8  (93) 

assuming  C  has  a  continuously  turning  tangent.  Altogether  (8?)  becomes 

${#*-!  «aprtf  + 1  b“P^  - 1  b^r  + 

c 

+  ^  +  •*/%  +  ^Apta)]6W  +  i^nanp80nJjds  (9W 

From  this  the  boundary  conditions  on  C  may  be  read  off.  They  aroJ 
prescribe  [tf*  -  |  e%Y  ♦  \  bW  -  |  bftP'  *  or  U8  (95) 

■  *S  +  ^np *  fc*V«>  «  w  <?6> 

"  ^nQnp  or  (97) 

For  the  alternative  theory  with  (83)  for  the  bending  strain  and 


t 


(P6;  (67)'  and  (865  for  equilibrium  equations,  the  boundary  conditions  aie: 

prescribe  [S^  -  |  +  bgCS^+t5^)^"^]^  or  Uq  (95)' 

"  At  *  +  I^Vp5  or  W  (96)' 

n  Tf^nn.  or  0  (97)' 

a  p  n 

FURTHER  APPROX IHAT IONS 
Snail  Rotation  About  the  Normal 

If  the  rotation  about  the  normal  can  be  neglected  compared  to  the  other 
two  rotations,  then  the  equations  can  be  simplified  further.  The  importance 
of  the  rotation  about  the  normal  is  not  entirely  established  at  the  present 
time  and  no  general  condition  under  which  it  can  be  neglected  is  known. 
Several  linear  theories  for  thin  shells  have  been  constructed  which  differ 
from  love’s  first  approximation  only  by  terms  in  the  bending  strain 
proportional  to  the  rotation  about  the  normal.  The  differences between  these 
theories  and  love's  are  tabulated  in  reference  11  where  the  general  validity 
of  these  theories  is  questioned.  That  the  rotation  about  the  normal  can 
sometimes  bo  neglected  is  evidenced  by  the  fact  that  these  theories  lead 
to  very  nearly  the  same  results  as  more  accurate  theories  in  some  specific 
applications.  See,  for  example,  reference  1$.  On  the  other  band  these 
theories  lead  to  erroneous  results  in  other  applications.  See  references 
16  and  17. 

For  those  cases  in  which  the  approximation  is  valid  the  strains  I  . 

°P 

and  given  by  equations  (71)  and  (82)  can  be  simplified  to  read 

V '  Ks*v  *  V  ♦  I  Vp  <»> 

•  Kn'V 


(99) 


The  corresponding  approximate  equilibrium  equations  (obtained  via  the  virtual 
work  principle)  are 


(100) 

^a  -  +  p  ■ 0 

(id) 

iff  -tf-o 

(102) 

and  the  boundary  conditions  are  to  prescribe: 

(If3*3  +  bgt6tY>f^]np  or  Ua 

(103) 

+  nJT*^  ♦  ^Of\np)  or  W 

(10W 

^Vp  or  K 

(105) 

When  these  equations  are  linearized  they  reduce,  essentially, 

to  those 

given  in  reference  10. 

The  Donnell -Mushtari-Vlasov  Approximation 

A  further  simplification  of  the  above  equations  is  possible  under  the 

assumptions  discussed  in  reference  13.  This  consists  in  neglecting  the 

term  containing  Uq  in  the  expression  for  with  the  following  results 

for  strains, 

V^V'VrVi 

(106) 

Kop  "  a^Q,p^p,8^  "  "W,op 

(107) 

ra  -  K  +  w,a 

(108) 

for  equilibrium  equations, 

(109) 

*?«  - v®  *  *  p  - 0 

(no) 

at  -V-o 

>p 


and  for  boundary  conditions  prescribe 


tf^n,  or  U 
P  « 


°r  W 


a^n  n.  or  0 
g  p  n 


(in) 

(112) 

(113) 

(11W 


Marguerre's  Shallow  Shell  Equations 

If  applied  to  a  shallow  shell  the  preceding  equations  can  be  further  > 
simplified  because  of  the  geometry.  Suppose  that  the  shell  is  nearly  flat 
and  parallel  to  the  t?  •  z  •  0  plane,  and  that  the  squares  of  the  slopes 
of  the  shell  with  respect  to  the  z  -  0  plane  may  be  neglected.  Then, 
approximately! 


•  -2 


op 


(115) 


Since  the  displacements  are  considered  small  compared  to  W  ,  the 
horizontal  displacements  and  the  vertical  displacements  B  are  given 
approximately  in  terms  of  and  W  by 


TJ  « tT  +  z  B 
a  a  ,a 

V*  B 

In  terms  of  (T  and  W  the  membrane  strain  E 

<K  u 


(116) 

(equation  (106))  beoomesi 


(117) 


The  strain  and  xa  are  as  before  (in  equations  (107)  and  (108)),  The 
conditions  of  equilibrium  in  the  horizontal  and  vertical  directions  are 


+  P 


(118) 
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•  « 


(119) 


-25- 


tfj?  -  Tf  =  0  (120) 

>P 

where  pc  is  the  horizontal  load  intensity  and  p  is  the  vertical  load 
intensity.  The  boundary  conditions  are  to  prescribe 


or  0* 
ts 

(121) 

+  or  w 

(122) 

or  0n 

(123) 

These  are  Marguerre's  shallow  shell  equations  in  tensor  form  (reference  1). 

CONCLUDINQ  REMARKS 

Several  nonlinear  theories  for  thin  shells  have  been  derived  in 
increasing  stages  of  approximation.  The  linearization  of  these  equations, 
which  is  more  or  less  obvious,  has  been  omitted  but  in  most  cases  the 
resulting  linear  equations  are  essentially  the  same  as  shell  equations 
already  given  in  the  literature.  In  all  cases  the  theories  are  first 
approximation  theories  in  the  sense  that  transverse  shear  and  normal  strains 
are  neglected. 

In  each  of  the  theories  derived  in  this  paper  the  equilibrium  equations 
and  strain-displacement  relations  are  related  by  a  principle  of  virtual  work 
and  hence  the  usual  variational  principles  may  be  formulated  and  proved. 
These  derivations  are  also  more  or  less  obvious  and  have  been  omitted. 

The  additional  manipulations  necessary  to  apply  the  equations  to 
stability  problems  has  not  been  given  either,  but  the  process  is  well  known 
and,  of  course,  different  manipulations  may  be  required  in  different  special 


cases. 


SHELL  EQUATIONS  IN  ORDINARY  NOTATION 
In  the  ordinaiy  notation  with  lines  of  curvature  for  coordinates  (as 
used  in  references  9  and  lli)  the  expressions  (71)  for  middle  aurface  strains 
and  (82)  for  bending  strains  are,  respectively 

1/ 1  3U2  1  3U1  1  3al  „  1  3a2  ..  . 

e12  ’  2(^3^  +  ^3^-^3^Ul  "^^U2  +^2) 

1  3?*!  1  9«1 

1  302  1  3c2 

*22  ‘  “2  3^  +  V2  351  ^1 


(A-l) 


(A-2) 


„ ifi. +  j, +  i/i  „ 

12  ^  3^  a2  3 ^  OjQg  3^  ™1  3^  ^2  2^J£  E™ 


\a2352 

The  transverse  shear  strains  are  given  by 

v  .  1. 3W  U1  v 

rl  5V5  * 


(A-3) 
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The  rotations  are  given  by 
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The  equilibrium  equations  (86),  (87),  and  (88)  now  read 
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The  boundary  conditions  on  an  edge  ■  constant  are  to  prescribe 
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The  boundary  conditions  on  an  edge  "  constant  are  the  same  as  these  with 
the  subscripts  1  and  2  interchanged  except  for  the  second  one  in  which  the 
sign  of  the  term  involving  0  should  be  changed. 

The  terms  in  the  preceding  equations  which  drop  out  when  rotations 
around  the  normal  are  neglected  hare  a  solid  underline.  The  terms  which 
drop  out  in  the  Donnell -Hushtari-Vlasov  approximation  are  those  with  either 
a  solid  underline  o.*  a  dashed  underline. 

The  bending  strains  of  the  alternative  theory  (equation  (83))  in  the 
present  notation  are 
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The  membrane  strains  are  the  same  as  those  given  by  (A-l).  The  force 
equilibrium  equations  are 
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The  moment  equilibrium  equations  are  the  same  as  (i-8)  and  (A-9). 

On  an  edge  ^  ■  constant  the  boundary  conditions  are  to  prescribe 
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On  an  edge  £g  -  constant  the  boundary  conditions  are  the  same  as  these  with 
the  subscripts  1  and  2  interchanged  except  in  the  second  line  the  term  with  0 
should  be  changed  in  sign. 
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